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ABSTRACT 

[OIII] A5007A emission lines of 16 intrac luster planetary nebulae candidates in the Coma cluster were detected 
with a Multi-Slit Imaging Spectroscopy (MSIS) technique using FOCAS on the Subaru telescope. The identi- 
fication of these faint emission sources as PNe is supported by (i) their point-like flux distribution in both space 
and wavelength, with tight limits on the continuum flux; (ii) the identification of the second [OIII]A4959 line in 
the only object at high enough velocity that this line too falls into the filter bandpass; (iii) emission line fluxes 
consistent with PNe at 100 Mpc distance, in the range 2.8 x 10"'^- L7 x 10"'^ ergs"' cm~^; and (iv) a narrow 
velocity distribution approximately centered on the systemic velocity of the Coma cluster Comparing with the 
velocities of galaxies in our field, we conclude that the great majority of these candidates would be intracluster 
PNe, free floating in the Coma cluster core. Their velocity dispersion is ^ 760km s"', and their mean velocity 
is lower than that of the galaxies. The velocity distribution suggests that the intracluster stellar population has 
different dynamics from the galaxies in the Coma cluster core. 

Subject headings: (ISM:) planetary nebulae: general; galaxies: cluster: general; galaxies: cluster: individual 
(Coma cluster); galaxies: evolution 



1. INTRODUCTION 

Cosmological simulations of structure formation predict 
that galaxies are dramatically modified by galaxy interactions 
during the assembly of galaxy clusters, losing a substantial 
fraction of their stellar mass which today must be in the form 
of intracluster stars (Murante et aJ. 2004). Observations now 
show that there is a substantial amount of intracluster stel- 
lar population, which is observed as diffuse intracluster light 
(ICL, Bernstein et al. 1995), or as individual stars, i.e., plan- 
etary nebulae (PNe; Arnaboldi et al. 2003, Feldmeier et al. 
2004) and red giant stars (DuiTell et al. 2002). The ICL rep- 
resents about 10% of the stellar mass in the Virgo cluster and 
as much as 50% in rich Abell clusters (Arnaboldi 2004). 

Intracluster planetary nebulae (ICPNe) are the best suited 
tracers for dynamical studies of the ICL because of their 
strong [OIII] 5007A emission, which allows both identifica- 
tion and radial velocity measurement. Measuring the pro- 
jected phase-space for the ICPNe constrains how and when 
this light originated (Napolitano et al. 2003). So far, ICPNe 
velocities have been measured only in the Virgo cluster at 15 
Mpc distance (Arnaboldi et al. 2004). Here we extend these 
studies to the Coma cluster (Abell 1656), one of the best stud- 
ied nearby galaxy clusters. Of these, it is the richest and the 
most compact one, providing a laboratory of prime impor- 
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tance for studying the effect of a dense environment on galaxy 
evolution. The ICL in Coma has been suspected since the 
early study by Zwicky in 195 1 and was confirmed by Thuan & 
Kormendy (1977) and Bernstein et al. (1995). The brightest 
PNe in the Coma cluster at distance of ~ 95Mpc (Bernstein 
et al. 1995) have fluxes of 2.2 x 10"'^ erg s"' cm"^ (Ciarduflo 
et al. 2002). To detect even the brightest of the Coma PNe, 
we must find a way to decrease substantially the noise from 
the night sky. 

We can achieve this with a technique that is similar to the 
approach used by Tran et al. (2004) and Martin & Sawicki 
(2004) to search for Lya emitting galaxies at very high red- 
shift. It combines a mask of parallel multiple slits with a nar- 
row band filter, centered around the [OIII]A5007A line at the 
redshift of the Coma cluster, to obtain spectra of all PNe that 
happen to lie behind the slits. Because the [OIII] emission 
lines from PNe are only a few km s"' wide, their entire flux 
still falls into a small number of pixels in the 2D-spectrum, 
determined by the slit width and seeing. On the contrary, the 
sky emission is dispersed in wavelength, allowing a large in- 
crease in signal-to-noise (S/N). The narrow band filter limits 
the length of the spectra on the CCD, so that more slits can 
be exposed. For brevity we will refer to this as the Multi Slit 
Imaging Spectroscopy (MSIS) technique. No conventional 
imaging technique can decrease the sky surface brightness in 
a similar way. In this Letter we show that with this technique 
we are indeed able to detect PNe in the Coma cluster, and 
briefly describe the first results obtained from one MSIS im- 
age in Coma. 

2. OBSERVATIONS 

The observations were carried out with the FOCAS spec- 
trograph at the 8.2m Subaru telescope on April 23, 2004. The 
spatial resolution of FOCAS is O'Mpix"', so the 6' diameter 
of the circular field of view (FOV) corresponds to 3600 pixels. 
We used grating 300B, which gives a measured dispersion of 
1.4-1. 5 A per pixel on the two FOCAS CCD chips. We used 
the N512 filter with FWHM of 60 A, centered at A, = 5 121 A, 
the wavelength of the redshifted [OIII] emission from a PN at 
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Fig. 1. — Multi-slit mask used in the observations, superposed on the 
observed field. North is to the left and East is up. The extended Ught at the 
NNW margin is from the halo of NGC 4874. The diameter of the FOCAS 
FOV, indicated by the circles, is 6' (about 165 kpc at Coma). The dispersion 
direction is West-East. Slitlets are arranged along a total of 70 long slits, but 
with short interruptions for mechanical stability of the mask. 



the mean velocity of the Coma cluster. The FWHM includes 
only ±1.6x the galaxy velocity dispersion in the Coma core. 
We will be able to detect the redshifted [OIII] 4959A emission 
only for the brighter PNe with the largest receding velocities, 
> 7400 kms"' , so that it can be seen near the blue edge of the 
filter. 

The light passing through the narrow band filter and a 0".6- 
wide long slit, and dispersed by the grism projects down 
to a spectrum of about 43 pixels. A mask was therefore 
constructed with uniform long slits spaced every 50 pixels, 
and interrupted only by short sections to ensure mechani- 
cal stability (see Fig. 1). The slitlet width on the mask 
was 0."6, corresponding to six 0."1 pixels. The area sur- 
veyed by this mask configuration is then about (3600 /50) slits 
x0".6 X 360" X (7r/4) = 12215 arcsec^, or 12% of the whole 
FOCAS FOV. 

With this mask, grism, and filter we took six 30 min ex- 
posures of a field centered at a =12:57:17, 5 =-1-28:09:35 
(B1950). This field is near the X-ray centroid of the Coma 
cluster and is essentially concident with the field observed 
by Bernstein et al. (1995). The measured seeing of the im- 
ages varied from 0."6 to 0."8. The CCD readout of the im- 
age was done with full spectral resolution and with a spa- 
tial binning of 2 pixels. Monochromatic, point-like emitters 
appear as elongated ellipses on the CCD, with width of ap- 
proximately 4 rebinned pixels and height of 5 pixels in the 
wavelength direction. Then the effective spectral resolution 
is ~ 7.3A, or 440kms"'. Data reduction was carried out in 
IRAF, with flux calibration using the spectrophotometric stan- 
dard star BDH-33d2642. 

2.1. How many PNe do we expect per slit configuration? 

We have estimated the average surface brightness of the dif- 
fuse light in our Coma core field from the results of Bernstein 
et al. (1995) to be = 24.7 mag arcsec"^. Assuming B-R = 
1.0, the average surface brightness in B in this field is 25.7. 
The total magnitude in the FOCAS circular field of 6' diame- 
ter is then Bjot = 13.18, which amounts to a total luminosity 
of Lb = 7.6 X 1O"*L0. To estimate the corresponding num- 



ber of PNe, we assume that the planetary nebulae luminosity 
function (PNLF) can be described by the formula of Ciardullo 
et al. (2002) which is a good fit to numerous observations in 
nearby galaxies. With the observations described below, we 
will be complete for PNe ^1.4 mag fainter than the bright 
cutoff of the PNLF. Using the luminosity-specific PN den- 
sity determined for an evolved population such as in M87, 
ai.o,B = 5.6 X 10"'' PN L^'q (Jacoby et al. 1990, Peletier et al. 
1990, Ciardullo et al. 1998), we expect 425 PNe associated 
with the diffuse light in this field, down to 1 mag fainter than 
the cutoff. 

The fraction of the FOCAS FOV surveyed by the mask used 
in our observations is 12% (see Section 3), giving about 50 
PNe located behind the mask slitlets. However, because the 
seeing FWHM is nearly equal to the slit width, of order half of 
these will have a significantly lower flux measured than their 
true flux; for FWHM=slitwidth, a point source at the center 
(edge) of the slit is dimmed by 24% (51%). Also, we lose a 
small fraction of the remaining PNe due to the limited filter 
bandwidth, if these PNe have the same velocity distribution as 
the Coma galaxies. Thus finally we expect to detect approxi- 
mately 20-30 PNe per mask. 

2.2. The influence of observational parameters on the 
detectability of PNe 

Here we show how the S/N of the detected PNe depends on 
the observational parameters, mainly seeing, slit-width, and 
spectral resolution. For a distance modulus of 34.9, a PN at 
the bright end of the PNLF has a flux in the [OIII] 5007A 
line of F5007.cutoff = 2.2 x 10""* erg sec"' cm"^. 0."6 resolution 
corresponds to 300 pc at this distance. 

For a sky surface brightness /zymag/arcsec^, the magni- 
tude of sky in a pixel a'Jxa'y my sky = - 2.51og((Tvf y), 
and the monochromatic flux from sky at 5500A in 
the pixel is \ogFsky,x = -0.4m,,,5i,. - 8.4, or Fsky.x = 4x 
10"^10"" '*^''o'^-o'yerg sec"' cm"^A"'. For a dispersion d and 
slit width 6", the true spectral resolution is dnue = dniO/ay), 
where k is a constant of order unity. In our observa- 
tional set-up, the CCD is rebinned 2x1, giving a spa- 
tial resolution ct^xct," = 0".2x0."l, d = 1.45A/pixel, 9" = 
0."6, and k = 5/6. Thus the flux from the sky sur- 
face brightness in the rebinned FOCAS pixel is fsky.pix = 
4x 10~'^10~^'^'^^ axCTydtrueSi'g sec"' cm"^. For our instrumen- 
tal set-up, we have measured on the flux-calibrated two- 
dimensional spectra fsky,pix = 7.7 x 10"'^erg sec"' cm"^, which 
gives an effective /Lty = 22.2at 5120A, inserting our pixel sizes 
and spectral resolution. 

The PN as a point-like, monochromatic source falls onto 
n 2± k' iKfp- /{Aaxdy) pixels, where in our observations the see- 
ing is = FWHM = 0".6 and k' = 1.4. We can thus com- 
pute the total signal-to-noise for a PN source. For a total 
integration time of 3 hrs, i..e. t^xp = 1.08 x 10'* sec, a total 
telescope area of Stei = 5.281 x 10^ cm^, and an overall ef- 
ficiency of telescope -i- spectrograph -i- airmass, e = 0.1, and 
hv = 3.61 x 10"'^erg = E^^^^Q at A = 5500A we obtain 

SNRpw = 7.0 X 10""-4'^'« (1) 

J f )-"'(J^\'(±]-"'(I^V" 

Vl.45A/pixy V0."6y \6ayJ \3hrs) 

for a PN fainter than the cutoff by Am magnitudes. Relative 
to the sky noise in one pixel, the S/N of a PN at the PNLF 
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Fig. 2. — Two-dimensional median-averaged spectra of emission ob- 
jects in the FOCAS field. Wavelength is along the vertical axis (507.7- 
514.15 nm) with a rebinned resolution of 1.5A per pixel; the true spec- 
tral resolution is 7.3A or 440kms"'. The horizontal direction is along 
the mask slitlets, with a rebinned resolution of 0."2 per pixel. The left 
four panels show PNe candidates from the brightest to one of the faintest 
in the field. The fluxes are (34.7,31.9,18.6,5.6) ADU, corresponding to 
(17, 16,9, 3) X 10"" ergs"' cm"~. The rightmost panel shows the spectrum of 
a background galaxy with continuum and strong absorption, probably blue- 
wards of Lyo, and a possible line emission, probably Lyo. 



cutoff (Am = 0) would be 32. Both values assume that the 
emission from the PN falls through the slit completely. De- 
tection of a source is normally considered to be secure if the 
S/N relative to the sky noise in one pixel is > 9. In our data 
this corresponds to SNRpat > 2, i.e., 95% probability of the 
source being real. If we use this criterion for secure detection, 
we can detect PNe about 1.4 magnitudes down the PNLF. 

Equation (1) shows that independent of all other parame- 
ters, good seeing is critical for these observations. Also the 
SNR/JAf increases like the inverse square root of the true spec- 
tral resolution. However, so that the light of a PN falls through 
the slitlet completely, we must have > (/>. Because large 
slitwidths degrade the spectral resolution, there exists a trade- 
off between SNRp^r and the number of PNe that can be ob- 
served throught the mask (previous subsection). Consider that 
the PNLF has a sharp cutoff and then it is nearly flat. There- 
fore, as long as the SNRp^ remains large enough that we can 
detect PNe one magnitude fainter than the bright cutoff, it is 
advantageous to use larger slit-widths. For with a larger slit- 
width, the additional PNe that come through the mask out- 
number the PNe that are lost at the faintest fluxes. Typically, 
to observe one slit configuration for 5-6.5 hours will work as 
long as the seeing (/> is not worse than 0."8. 

3. RESULTS AND DISCUSSION: PNE IN THE COMA CLUSTER 

A total of 16 spectrally unresolved, point-like sources with 
undetected continuum were detected on the two CCD chips. 
Their emission line fluxes were in the range 2.8 x 10"'''- 
1.7 X 10"'^ ergs"' cm"^. Objects with a flux of less than 
2.5 X 10"''' ergs"' cm"^ (5 ADU) were considered as not real. 
Figure 2 shows two-dimensional spectra of the two bright- 
est objects, an intermediate flux source, and a low S/N emis- 
sion. These are all PN candidates. Fig. 2 also shows a source 
with a possible emission line and a continuum redwards of the 
line, but no continuum on the blue side, most likely a Lya- 
absorbed background galaxy. 

We now consider the evidence that the unresolved emission 
sources are in fact planetary nebulae. Figure 3 shows their 
1-d spectra, obtained by summing the respective columns in 
the 2-d-spectra. The intermediate flux source in Fig. 2 is at 
sufficiently large recession velocity that, if the observed line 
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Fig. 3. — One-dimensional spectra of the PNe candidates and of the proba- 
ble Lya galaxy shown in Fig. 2 (same order from bottom to top as in Fig. 2). 
The top panel shows the arclamp image with the two emission lines visible 
in our spectral range, observed with the same set-up. The widths of the ar- 
clamp lines are very similar to the widths of the brighter emission lines of the 
PNe candidates, showing that these are uni'esolved in wavelength. The object 
in the third panel from below is sufficiently redshifted that both lines of the 
oxygen doublet [OIII]A4959A, A5007A fall in the filter bandpass. The S/N 
of the emissions from the four PN candidates shown are (20, 18, 1 1, 3.3) with 
respect to the noise per pixel of 1.7 ADU or 8.5 X 10"-" ergs"' cm"^. 



is [OIII]A5007A, the second A4959A line also falls into the 
narrow band filter In fact, its 1-d-spectrum in Fig. 3 does 
show both lines of the doublet: the ratio of the equivalent 
widths of both lines is consistent with the theoretical value 
of 1:3. 

The brightest object has measured counts of 34.7 ADU in 
the median image, corresponding to a total line flux of 1 .7 x 
10""^ erg s"' cm"2, and is detected at S/N= 20. The faintest PN 
candidate has a S/N= 3.3, a line flux of 5.6 ADU, correspond- 
ing to 2.8 X 10"'''ergs"' cm"^, or a total of 68 [OIII]A5007A 
photons in three hours, about one photon every three minutes. 
The flux from a PN in Coma at the bright cutoff of the PNLF 
is 2.2 X 10"'^ergs"' cm"^ (Section 2)^. Note that our emission 
sources are generally not centered in the slitlets, and even if 
they are, we lose 24% (38%) of the flux, given the slitlet width 
of 0."6 and a seeing FWHM of 0."6 (0."8). Thus the bright- 
est fluxes we measure are consistent with PNe in Coma at the 
cutoff of the PNLF 

Our candidates are also undetected in the continuum. 
The l(T-upper limit on the pixel flux for the extracted 
spectra in the range 5070A-5140A, excluding the line, 
is 3.4 ADU, corresponding to a continuum flux limit of 

1.6 X 10"^" ergs"' cm"^ A . Our brightest object thus has 
EW>110A. The V-band continuum flux of an 08 star at the 
distance of Coma is 4.3 x 10"^' erg s"' cm"^A - with our ob- 
serving technique we would see the continuum from a few O 
stars in the Coma cluster! This rules out compact HII regions 
such as or brighter than those observed in Virgo (Gerhard et al. 
2002). Similarly, background galaxies would often have a 

^ In the [OIII] PNe magnitudes defined by Jacoby (1989) these two objects 
have mjoo? = 30.7 and msoov = 32.5, and the bright cut-off of the PNLF in 
Coma is at mjoov cutoff = -4.5 + 34.9 = 30.4. 
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Fig. 4. — Velocity histogram for all 16 PNe candidates in the mask. Mean 
and median velocities are v = 6100kms"' and 

^'med — 6224 km s , respec- 
tively, and the velocity dispersion is cr = 760 km s ^ and cTmed 

= 562kms-' 

with respect to the median. Overplotted is the measured filter transmission 
curve in our observing set-up (full red line), using the continuum light from 
the spectrophotometric standard star BD+33d2642, and (dashed curve) this 
function times a Gaussian with vcoma = 6853 km s"' and crcoma = 1082 km s"' 
as for the galaxies (Colless & Dunn 1996). 



much larger continuum flux than our Hmit; for example, see 
the background object shown in Fig. 2. For our brightest can- 
didates, the large equivalent widths are larger than those of 
known [Oil] emitters at z = 0.37 (Hogg et al 1998). The line 
fluxes we measure are 10 - 100 times fainter, and the surface 
number density is 10 - 100 times higher, than for currently 
studied Lya emitters at z ^ 3 (Fynbo et al. 2003). Moreover, 
in the flux decade brighter than our brightest line flux, we see 
no object without measurable continuum, and only one object 
with blue continuum, for which EW=50A. 

Figure 4 shows the velocity histogram for our PNe candi- 
dates. They have velocities between 4800-7560kms~'. With 
our instrumental setup, we detect wavelengths in the range 
5071-5135A (FWHM), or velocities in 3830 -7680 km s"'. 
The velocity histogram is clearly not uniform in the filter 
bandpass (Fig. 4); an unclustered population of background 
emission line galaxies would appear as such a nearly uniform 
distribution, because only the faintest objects near the filter 
edges would fall out of the sample. 

In summary, we have the following evidence that our un- 
resolved emission sources are PNe in the Coma cluster: (i) 
They are unresolved both spatially and in wavelength, (ii) 
They have no detectable continuum, to a limit of 1.6 x 



10 ^"ergs 'cm The brightest object thus has EW> 

11 OA. (iii) We have seen both lines of the [OIII] doublet in 
the only source at sufficiently large recession velocity that 
also A4959A is redshifted into the wavelength range probed, 
(iv) The emission fluxes of the brightest objects are consis- 
tent with those of the brightest PNe in a population at dis- 
tance 100 Mpc. (v) The number density of our candidates 
is consistent with that expected from the measured surface 
brightness of the ICL in our Coma field, (vi) The distribu- 
tion of recession velocities is centered around the Coma clus- 
ter and is inconsistent with a population of background ob- 
jects uniformly distributed in velocity. However, we cannot 
rule out that a population clustered in velocity of so far un- 
detected low-luminosity background emission galaxies could 
contribute to our sample, if these exist in sufficient numbers 
at z - 3 (Lya) or z - 0.37 ([Oil]). 

From Bernstein et al. (1995), approximately half of the 
total luminosity in our field is in the ICL component of the 
Coma cluster, the other half is in the brighter galaxies, while 
faint galaxies do not contribute significantly. In the PNe, we 
expect a larger diffuse light fraction: PNe near bright galaxies 
will be harder to see because of the increased background. We 
have superposed our PNe candidates on the image of the field 
and compared their velocities with those of the brighter galax- 
ies in the field. None of our candidates is close in space and 
velocity to any of these, but two PNe have measured veloci- 
ties and positions that could be consistent with the cD enve- 
lope of NGC 4874, which is just outside our field, ^ 120 kpc 
from the field center NGC 4874 itself has a radial velocity 
of 7224kms"' and a velocity dispersion of 280kms"' (Smith 
et al. 2000). Most of our PNe, which have v = 6100kms"' 
and a = 760km s"', are therefore not associated with the halo 
of NGC 4874. 

Thus, the great majority of the detected PNe candidates 
would be intracluster PNe (ICPNe) in the Coma cluster. It is 
interesting that the measured mean and dispersion velocities 
for these ICPNe are smaller than those measured for the Coma 
cluster core centered around NGC 4874, vcoma = 6853 kms"', 
velocity dispersion crcoma = 1082kms"' (Colless & Dunn, 
1996). A KS test shows that the ICPN velocity distribution 
has probability 12% to be drawn from the expected veloc- 
ity distribution of galaxies in our filter (Gaussian times filter 
band transmission). This suggests that the dynamics of the 
intracluster population in the Coma core may differ from that 
of the galaxies, and calls for further investigation. 

We acknowledge financial support by the Swiss SNF and 
by INAF. We are grateful to the on-site Subaru staff for their 
support. 



REFERENCES 



Amaboldi. M., 2004, in lAU Symp. 217, Recycling intergalactic and 

interstellar matter, eds. RA. Due , J. Brian, E. Brinks, ASP, 54 
Arnaboldi, M. et al. 2003, AJ, 125, 514 
Arnaboldi, M. et al. 2004, ApJ, 614, L33 
Bernstein, G.M., et al. 1995, AJ, 110, 1507 
Ciardullo, R. et al. 1998. ApJ, 492, 62 
Ciardullo, R., et. al. 2002, ApJ, 577, 31 
Colless, M., & Dunn, A.M. 1996, ApJ, 458, 435 

DuiTell, P. R., Ciardullo, R., Feldmeier, J. J., Jacoby, G. H., & Sigurdsson, S. 

2002, ApJ, 570, 119 
Feldmeier, J. J., Ciardullo, R., Jacoby, G. H., & Durrell. P R. 2004, ApJ, 615, 

196 

Gerhard, O., Amaboldi, M., Freeman, K.C., Okamura, S. 2002, ApJ, 580, 
L121 



Hogg, D.W., Cohen, J.G., Blandford, R., & Pahi'e, M.A. 1998, ApJ, 504, 622 
Jacoby, G.H. 1989, ApJ, 339, 39 

Jacoby, G.H., Ciardullo, R., Ford, H.C. 1990, ApJ, 339, 39 
Martin, C.L., & Sawicki, M. 2004, ApJ, 603, 414 
Murante, G. et al. 2004, ApJ, 607, L83 
Napolitano, N.R. et al. 2003, ApJ, 594, 172 

Peletier, R.R, Davies, R.L., Ilhngworth, G.D., Davis, L.E., Cawson, M. 1990, 
AJ, 100, 1091 

Smith, R.J., Lucey, J.R., Hudson, M.J., Schlegel, D.J., Davies, R.L. 2000, 

MNRAS, 313, 469 
Thuan, T.H., Kormendy, J. 1977, PASP 89, 466 

Tran K.-V., Lilly. S.J., Crampton, D., Brodwin, M. 2004, ApJ, 612, L89 



Fynbo, J.PU., et al. 2003, A&A, 407, 147 



